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Coincident observations with gravitational wave (GW) detectors and other astronomical instru- 
ments are in the focus of the experiments with the network of LIGO, Virgo and GEO detectors. 
They will become a necessary part of the future GW astronomy as the next generation of advanced 
detectors comes online. The success of such joint observations directly depends on the source lo- 
calization capabilities of the GW detectors. In this paper we present studies of the sky localization 
of transient GW sources with the future advanced detector networks and describe their fundamen- 
tal properties. By reconstructing sky coordinates of ad hoc signals injected into simulated detector 
noise we study the accuracy of the source localization and its dependence on the strength of injected 
signals, waveforms and network configurations. 

PACS numbers: 04.80.Nn, 07.05. Kf, 95.55.Ym, 04.30. Db 



I. INTRODUCTION 

There has been a significant sensitivity improvement of 
the gravitational wave detectors since the Laser Interfer- 
ometer Gravitational Wave Observatory (LIGO) pQ and 
Virgo observatory [2 started their operation. In 2007 
LIGO and Virgo completed the two year run at sensitiv- 
ity that allows detection of a merger of two neutron stars 
(NS-NS) as far as - 30 Mpc away j3H]. In the most re- 
cent run (May 2009 - October 2010) the binary neutron 
star horizon distance has been increased to ~ 40 Mpc. 
However, even at this impressive sensitivity, the antic- 
ipated detection rate with the initial LIGO and Virgo 
detectors is quite low. A detection may be possible in 
the case of a rare astrophysical transient event such as a 
supernova explosion in our Galaxy or a nearby merger of 
binary neutron stars. The signal is likely to be weak and 
it will be difficult to prove its astrophysical origin un- 
less it is confirmed with a coincident observation of the 
electromagnetic or neutrino counterpart. For this rea- 
son the LIGO and Virgo collaborations are conducting a 
wide range of joint observations [5] with other astrophys- 
ical experiments including radio [6j [7] , optical and x-ray 
telescopes [8UTT] , and neutrino detectors [I2j [13] . 

A more robust detection of gravitational waves from 
astrophysical sources is anticipated in the next five years 
as Advanced LIGO and Advanced Virgo come online. 
Numerous GW signals, expected to be observed by ad- 
vanced detectors (likely ~ 40 NS-NS events per year [T4]), 
will begin our exploration of the gravitational- wave sky 
and start the era of the gravitational wave astronomy. 
Along with the advanced GW detectors, a new generation 
of optical telescopes will come online [T5HT7] , which will 



enable a wide and deep survey of the electromagnetic sky. 
Joint observations with the advanced gravitational wave 
detectors and electromagnetic instruments will not only 
increase the confidence of detection but also bring funda- 
mentally new information about the GW sources. They 
will reveal the physics and dynamic of sources, provide 
the identification of host galaxies and the associated red- 
shifts, and in some cases determine luminosity distance 
to the source. 

One of the major challenges for such joint observations 
is to establish unambiguous association between a gravi- 
tational wave signal and a possible electromagnetic coun- 
terpart. It greatly depends on the ability of the GW net- 
works to reconstruct sky coordinates of a detected GW 
source. Given an accurate sky location, a correspond- 
ing electromagnetic transient may be identified in a list 
of events obtained with the all-sky telescope surveys, or 
the EM instruments can be guided to take images of a 
small area in the sky. In the second case, it is important 
that the sky localization is performed by GW detectors 
in real time with low latency. The efficiency of the GW- 
EM association and the choice of a partner telescope is 
affected by the sky localization error which should be 
well within the instrument's field of view (typically less 
than few square degrees). Moreover, exploring smaller 
area in the sky will decrease the probability of the false 
association. 

The problem of the source localization with networks 
of GW detectors is in the focus of research in the gravi- 
tational wave data analysis. There are several analytical 
studies [T8H2T] of this problem considering geometrical 
reconstruction of source coordinates based on the tri- 
angulation, which requires a measurement of the arrival 
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time of a GW signal at different detectors. However, the 
accurate timing of the GW signal is intimately related to 
the reconstruction of the signal waveforms. Due to the 
different detector sensitivities to the GW polarizations, 
the waveforms recorded by individual detectors may be 
different and they may not have a common timing refer- 
ence (like a signal peak time) for a direct measurement of 
the differences in the arrival time. Therefore, the problem 
of the source localization is better addressed in the frame- 
work of the coherent network analysis [22-24J, which re- 
constructs the waveforms and the sky coordinates simul- 
taneously. By using both these methods (triangulation 
and coherent network analysis), several practical source 
localization algorithms [25 -27 have been recently devel- 
oped and used during the LIGO and Virgo data taking 
runs in 2009-2010. 

There have been a number of studies addressing ben- 
efits of individual detectors [28j [29] and various detector 
networks [30] [31] . In this paper we present a simulation 
study of the source localization and the reconstruction of 
GW waveforms with the networks of advanced detectors. 
The study is performed with a coherent network method, 
called coherent WaveBurst [25 (cWB), based on the like- 
lihood analysis . In cWB the data from all detectors in 
the network is processed simultaneously in order to re- 
construct a common GW signal which is consistent with 
the recorded detector responses. The consistency is mea- 
sured by the likelihood ratio, which is a function of the 
source parameters (waveforms and sky location). The 
most probable source parameters are obtained by maxi- 
mizing the likelihood ratio over the signal waveforms and 
sky coordinates. The method performs reconstruction of 
unmodeled burst signals (arbitrary waveforms) and sig- 
nals with a certain polarization state: elliptical, linear 
and circular. 

The paper is organized as follows. Possible networks 
of advanced detectors and their fundamental properties 
are descussed in section [ill In sections IIIII we describe 
the reconstruction algorithm. The simulation framework 



for this study is presented in section |TVj The results are 
reported in section |V| In sections [VT| and \VU\ we describe 
main factors limiting the source reconstruction and dis- 
cuss the results. 



II. DETECTOR NETWORKS 

In 2001-2010 the LIGO Scientific Collaboration (LSC) 
and the Virgo Collaboration operated a network of in- 
terferometric gravitational- wave detectors which are the 
most sensitive instruments from the first generation of 
the GW interferometers (1G). They consist of power- 
recycled Michelson interferometers with kilometer-scale 
Fabry-Perot arms designed to detect gravitational waves 
with frequencies between tens of Hz and several kHz. 
The two LIGO observatories [T] are in Hanford, Wash- 
ington (4 km and 2 km detectors) and in Livingston, 
Louisiana (4 km detector), and the 3-km Virgo detec- 
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TABLE I: Geographical locations and orientations of the 2G 
detectors. The orientation of the detector arms is defined by 
the rotation angle (counterclockwise) with respect to the local 
coordinate frame with axises due North and East. 



tor [2] is located in Cascina, Italy. Other gravitational 
waves interferometers are the 300 m detector TAMA [32] 
in Mitaka, Japan, and the 600 m detector GEO600 [33] 
in Hannover, Germany. Currently all 1G interferometers 
are decommissioned, except Virgo and GEO600, which 
continue to take data. 

The second-generation GW detectors (2G) are cur- 
rently under construction. They include the advanced 
LIGO detectors [34], and the advanced Virgo detector 
(V) [35] which will have by an order of magnitude better 
sensitivity than the 1G detectors. All advanced LIGO 
detectors have 4-km long arms, with one detector in Liv- 
ingston (L) and two identical co-aligned detectors in Han- 
ford (H and H). Also there are plans to build the Large 
Cryogenic Gravitational Telescope (LCGT) [36] [37] in 
Japan (the J detector) and possibly move the LIGO H 
detector to a site in Australia [30] [38] (the A detector). 
Figure [I] shows the design sensitivity for the listed 2G 
detectors. These, hopefully all five interferometers, com- 
pose the most advanced GW detector network which will 
be in operation after 2015. 




Frequency (Hz) 



FIG. 1: Amplitude spectral density of the design noise for the 
second generation detectors. 

The network performance greatly depends on the num- 
ber of detectors in the network, their location and orien- 
tation of the detector arms. Table [I] shows the geograph- 
ical coordinates of the instruments and the orientation 
of the detector arms used in this study. For the Aus- 
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tralian instrument the orientation of the detector arms is 
not yet decided, therefore we consider two possible con- 
figurations: A - arms are due north and east, and A - 
arms are rotated counterclockwise by 45° with respect to 
A. Depending what instruments are constructed we con- 
sider several network configurations, as listed in Table \H\ 



A. Network sensitivity 

The sensitivity of the network of K detectors is fully 
characterized by its noise-scaled antenna pattern vectors 
f + and f x : 



f +(x) 
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(2.1) 



where o\ is the variance of the noise and (Fk+,Fk X ) are 
the antenna patterns of individual detectors. In general, 
the <Jk is a function of frequency. For a discrete data 
with the Nyquist frequency /n in the Fourier or wavelet 
domain, the Sk = v\j ' Jn is the estimator of the single- 
sided power spectral density of the detector noise. The 
power spectral density S^et and the variance <j^ et of the 
network noise are defined by the following equation: 



K 



$net — 



^net 

In 



= US 



(2.2) 



k = l 



Note, that both the S^et and cr^et decrease as more de- 
tectors are added to the network. A network of K equally 
sensitive detectors has by a factor \/K lower noise am- 
plitude than the individual detectors. 

The antenna patterns depend upon the source coordi- 
nates (0, (j)) and the polarization angle which defines 
the wave frame of the incoming GW signal h = h x ). 
It is convenient to define vectors f + , f x and h in the dom- 
inant polarization wave frame [24] where (f + -f x ) = and 
|f+| > |f x |. The vectors f + and f x define a vector of the 
noise-scaled detector responses to the wave h 



(2.3) 



The inner product (£h|£h) calculated over the sampled 
detector responses gives the estimator of the network 
signal-to- noise ratio (SNR): 



(2.4) 



In general, the inner product of two network vectors a 
and b is defined as 



(a|b) = £>[i]-b[i]) 



(2.5) 



where the sum is taken over the data samples i containing 
the event. 

The norms of the antenna pattern vectors |f_|_| and |f x | 
characterize the network sensitivity to the GW polariza- 
tions. To illustrate this and other network properties, 



we assume below that in the signal frequency band the 
vectors f + and f x do not vary much. In this case 



Pnet « V|f+| 2 (^+IM + |fx| 2 (^x|M (2.6) 

where the inner product (h|h) = (/i + |/i + ) + (h x \h x ) de- 
termines the root-sum-square amplitude of the GW po- 
larizations: 
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As it follows from Equation |2.6[ the network alignment 
factor [25] 



|f+i 



(2.8) 



characterizes the relative network sensitivity to the two 
GW polarizations. It determines the ratio of the SNRs 
from each GW component, assuming that in average 
their sum-square energies are the same: (h+\h+) = 
(h x \h x ). Closely aligned networks (like LHH) have poor 
sensitivity to the second polarization (a <C 1) making 
difficult the reconstruction of the full GW signal. 

The overall network sensitivity is characterized by the 
effective power spectral density of the network noise 
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which depends on the sky coordinates. It determines the 
average network SNR for a population of GW signals 
with the average amplitude h rss /V2 per polarization: 



h r 



Pnet 



(2.10) 



It is convenient to factorize the sky-dependent part of 

A^net as 



A^net — ^ 2 Snet j 



(2.11) 



where T is the network antenna factor distributed be- 
tween (low sensitivity) and 1 (high sensitivity). Fig- 
ure [2] shows the antenna and the alignment factors for 
different networks as a function of the latitude and lon- 
gitude of the source (skymaps). Since these network pa- 
rameters are noise dependent, the skymaps are calculated 
at the frequency 100 Hz where the advanced detector sen- 
sitivities are about the same. The T distribution shows 
how uniform is the network response across the sky. The 
values of a close to unity indicate the same sensitivity 
to the two GW components. Respectively, the values of 
a close to zero indicate that the second GW component 
is not measurable for a weak GW signal. As Figure [2] 
shows, several non-aligned detectors (preferably five) are 
required for elimination of these blind spots in the sky. 

One of the main characteristics of a detector network 
is its search volume. Given an isotropic distribution of 
transient sources with the root-square-sum amplitude h Q 
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FIG. 2: T/ie distributions of T (left plots) and a (right plots) at the frequency of 100 Hz as a function of latitude (0) and 
longitude (<j>) for the networks LEEV, LEVA, LEVA, LEV J, LEVA J (from top to bottom). 

at the fiducial distance r G , the search volume is defined as [39] 
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300 Hz 
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1.63 


1.15 


1.80 



TABLE II: Expected difference in detection rates with respect 
to the LHV network. 



where e is the detection efficiency. Assum ing the same 
SNR thresholds p ne t(h) (see Equation 2.10), the V ne t can 
be calculated with respect to the volume Vo of the refer- 
ence network 
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where < N^t > and < 7V o/ > are the averages over 
the sky. The ratio V net /Vo is quite independent on the 
search algorithm and the GW source model. Table [TT] 
shows the volume (and detection rate) ratios calculated 
with respect to the LHV network (Vo = Vlhv)- As more 
detectors are added to the network the detection rates 
increase. Being beneficial this increase, however, is not 
critical for the first direct observation of gravitational 
waves. More important, the networks with more detec- 
tors are likely to be less affected by non- Gaussian and 
non- stationary noise than the networks with fewer de- 
tectors. They are expected to have lower false alarm 
rates and higher detection confidence for the same p net 
threshold. Also, re-location of the LIGO H detector to 
Australia does not affect the detection rates for both A 
and A configurations. But, as shown in Figure [5J the 
LHVA detector configuration would be more preferable 
for the reconstruction of the GW polarizations than the 
LHHV and LHAV networks. 



III. RECONSTRUCTION ALGORITHM 

A. Coherent network analysis 

One possible approach to the coherent network analy- 
sis is based on the Neyman-Pearson criterion which de- 
fines the likelihood ratio 



A(x,fi) = 



p(x\0) 



(3.1) 



where x is the network data vector, the p(x|0) is the 
joint probability that the data is only instrumental noise 
and p(x|h) is the joint probability that a GW signal h is 
present in the data. In general the likelihood ratio is a 
functional which depends upon the source parameters Q. 
One generalization of the Neyman-Pearson criterion is to 
maximize A(x, Q) over Q. The obtained maximum like- 
lihood ratio statistic reaches its maximum for the best 
match of the corresponding waveform to the data. If the 
source model allows calculation of the GW waveforms 
as a function of a small number of source parameters 



(for example, for binary black holes), then a template 
bank can be generated. In this case the variation is per- 
formed over the template bank and the likelihood ap- 
proach is equivalent to a matched filter. The cWB algo- 
rithm searches for unmodeled burst signals. In contrast 
to the binary black hole sources, where the number of 
parameters is relatively small, the parameters character- 
izing the unmodeled bursts are essentially the signal am- 
plitudes themselves at each instance of time. It is not 
possible to generate a template bank for such a large pa- 
rameter space. Instead, the best matching waveform is 
found by variation of A over unknown GW waveforms h. 



B. GW waveforms 

For stationary Gaussian noise the coherent WaveBurst 
algorithm defines the likelihood C as twice the logarithm 
of the likelihood ratio A 



£[h] = 2(w|60 - , 



(3.2) 



where the vector w represents whitened data from K 
detectors with uncorrelated noise 
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Xi[z,Ti] X k [i,T k ] 



X K [i,T K ] 



(3.3) 



The sampled detector amplitudes Xk[i,Tk] take into ac- 
count the time-of- flight delays ?>, which in turn depend 
upon the source coordinates 6 and <p. The solutions for 
the GW waveforms h, defined in the dominant polariza- 
tion frame, are foun d by the variation of the likelihood 
functional (Eq. pT2])): 

(3.4) 
(3.5) 



tf + = (w.f + )/|f + | 2 
H x =(wf x )/|f x | 2 



The maximum likelihood ratio statistic is calculated by 
substituting the solutions into C[h]. The result can be 
written as 



£max = ^2 wPwT ' 



(3.6) 



where the matrix P is the projection constructed from 
the components of the unit vectors e + and e x along the 
directions of the f + and f x respectively: 



Pfim — 6_|_ n 6_|_ r; 



(3.7) 



The kernel of the projection P is the signal plane defined 
by these two vectors. The null space of the projection 
P defines the reconstructed detector noise n rec which is 
referred to as the null stream. 

The projection matrix is invariant with respect to the 
rotation in the signal plane where any two orthogonal 
unit vectors can be used for construction of the P nm . 
Therefore one can select unit vectors u and v such that 
w • v = and then 

(3.8) 
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The unit vector u defines the vector 
£ = (w • u)u 



(3-9) 



whose components are the standard likelihood estimators 
of the noise-scaled detector responses 

C. Source coordinates 



The maximum likelihood ratio statistic L max is a func- 
tion of the sky coordinates and <j). If no information 
regarding the source coordinates is available then the 
variation over the sky should be also performed. It is ex- 
pected that L max takes maximum close to a true source 
location, however, it is not necessarily the optimal statis- 
tic for the coordinate reconstruction. The coherent part 
of the likelihood quadratic form 



n 7^ m 



(3.10) 



has a strong dependence on the time delays between the 
detectors and therefore the coherent energy E c is ex- 
pected to be a better statistics for the source localization. 
On the other hand the E c is a biased estimator: for an 
arbitrary GW signal it may take maximum value away 
from the true source location. To minimize the bias, the 
sky statistic is constructed in the following way 



-^sky — 



E(E - L n 



\Ec\) 



(3.11) 



where E = (w|w) is the total normalized energy of the 
signal and noise. This statistic penalizes the sky locations 
with low values of E c and large values of the null energy 
E — £max, and it reduces to the maximum likelihood 
statistic L max /E when the ratio (E — L max )/E c is small. 
The L s k y is used to rank different sky locations and calcu- 
late the probability distribution of the estimated source 
coordinates in the sky. 



D. Model-dependent constraints 

The likelihood method offers a convenient framework 
for introduction of constraints arising from the source 
models. Unlike for template searches where accurate 
waveforms are required, in principle, any useful informa- 
tion about sources can be used to constraint the likeli- 
hood functional. This allows customization of the generic 
burst algorithms in order to search for specific, but not 
very well modeled sources. One obvious class of con- 
straints is related to the different polarization states of 
the GW signals. For example, some of the core collapse 
models predict waveforms with a linear polarization [41] 
or random polarization [42]. Merging binary neutron 
stars or black holes are expected to produce elliptically 
polarized gravitational wave signals [43]. Also the neu- 
tron star mergers can be the source of the short GRB 



signals [44] where relativistic jets are emitted along the 
rotation axis of the binary system and in this case the as- 
sociated gravitational waves should have the circular po- 
larization. The cWB algorithm allows searches with sev- 
eral types of the polarization constraints: circular, linear, 
elliptical and random (or unmodeled search). All these 
searches are used in the study to estimate possible im- 
provement of the source localization if the reconstruction 
is constrained by the source model. 



IV. SIMULATIONS 

A. Injected signals 

Several types of "ad hoc" waveforms were used to study 
the performance of the detector networks for different 
signal frequencies and polarization states. They were in- 
jected into the simulated detector data streams in a wide 
range of signal-to-noise ratios with the coordinates uni- 
formly distributed in the sky. The Gaussian detector 
noise was simulated with the amplitude spectral density 
presented in Figure [I] The injected signals were band- 
limited white- noise waveforms (WNB) with the random 
polarizations and sine-Gaussian waveforms (SG) with the 
linear and circular polarizations. The WNB waveforms 
were bursts of white Gaussian noise in a frequency band 
(A 5 $2) which have a Gaussia n ti me profile with the stan- 
dard deviation r (see Table |TTT| ) . The random polariza- 
tion waveforms h+ and h x were selected to have the same 
square-sum energy: = (h x \h x ) 



Waveform 


r(s) 


fi (Hz) 


h (Hz) 


WNB 


0.1 


250 


350 


WNB 


0.1 


1000 


2000 



TABLE III: Simulated White Noise Bursts. 
The SG waveforms were simulated as follows 

h+{t) = h o sm(2irtf )exp(-t 2 /T 2 ) , 
h x (t) = fi! cos(27rf/ ) exp(-t 2 /r 2 ) , 



(4.1) 
(4.2) 



where fo is the waveform central frequency, h Q and hi are 
the waveform amplitudes, and r is related to the wave- 
form quality factor Q = y^irfor (see Table |lV| . The 
amplitude parameters were hi = for linear polarization 
and hi = h Q for circular polarization. During the anal- 
ysis the amplitude of injected events varied to simulate 
events with different signal-to- noise ratios. 



B. Error regions 

The injected signals are used for estimation of the accu- 
racy of the coordinate reconstruction. For each injected 
event the L s k y skymaps is calculated with the angular 
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Waveform fo (Hz) Q Polarization 
SGQ3 235/1053 3 Linear 
SGQ9 235/1053 9 Linear 
SGCQ9 1 235/1053 1 9 | Circular 



TABLE IV: Simulated sine- Gaussian waveforms with quality 
factors Q=3 and Q=9, low (235Hz) and high (1053Hz) fre- 
quencies, and two polarization types - linear and circular. 

resolution of dVt = 0.4 x 0.4 square degrees: ~ 200, 000 
sky locations (pixels) total. Figure [3] shows an example 
of such a skymap for one of the SGQ9 (235Hz) injec- 
tions. For such events it is typical to see a pattern of 




-180 -60 60 
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FIG. 3: Example of the likelihood sky map L s k y for an injected 
signal at = —30° and 4> = 144°: L s k y as a function of and 
cj) (top), L s ky distribution around the reconstructed location 
(bottom). 

fringes with large value of L s k y corresponding to a good 
match between responses due to a common GW signal 
reconstructed in different detectors. Such sky points are 
the most probable as the source location. Depending on 
many factors, such as the signal strength, waveform mor- 
phology, etc, the L s k y statistic can be well localized in a 
single small cluster in the sky or distributed over a large 
area which can be also split into several disjoint clusters. 
This type of ambiguity is typical for the least constrained 
unmodeled search and networks with only three spatially 
separated detectors. 

To characterize the accuracy of the coordinate recon- 



struction for a single injection we define an error region: 
total area of all pixels in the sky which satisfy the condi- 
tion L sky (#,(/>) > L sky (6i,(j)i), where (0*,0i) is the injec- 
tion sky location. Given a population of injected signals 
uniform in the sky, the 50 CL and 90 CL error regions, 
containing 50% and 90% of injections respectively, can 
be calculated. The median error angle is defined as the 
square root of the 50% error area. 

The L s k y skymap can be also converted into the prob- 
ability skymap which is normalized to unity if integrated 
over the entire sky. In this case the 50 CL and 90 CL 
error regions are represented by the most probable pixels 
with the cumulative probability of 50% and 90% respec- 
tively. Such probability skymaps are not relevant for the 
simulation studies we perform, but they are important 
for the analysis of real data. 



V. RESULTS 

A. Coordinate Reconstruction 

The accuracy of the coordinate reconstruction strongly 
depends on the strength of detected signals which can be 
conveniently characterized by the average (per detector) 
signal-to-noise ratio 

Pdet = Pnet/V^ . (5.1) 

For example, Figure [4] shows the dependence of the me- 
dian error angle 0^50% on pdet for all injected signals, 
which is well approximated by a function 

a 50% =A + B (—) + C (—) 2 . (5.2) 

\PdetJ \PdetJ 

The parameter A is the median error angle for events 
with very large SNR and A+B+C is the median error 
angle for events with pdet = 10. Figure [4] also shows a 
dependence of the coordinate resolution on the number 
of detector sites in the network. There is a significant 
improvement of the resolution when more sites are added 
to the network. This is particularly noticeable at low 
SNR, which is very important because the anticipated 
GW signals are likely to be weak. 

Because of several limiting factors (see section |VI[) 
the reconstruction is not uniform in the sky. Figure [5] 
shows the distribution of the median error angle across 
the sky for different network configurations. There is a 
dramatic improvement of the coordinate reconstruction 
for the LHVA, ALVJ and LHVAJ networks. However for 
the 4-site networks there remain areas where the source 
localization is poor. Figure [6] compares the pointing ca- 
pabilities of the network consisting of three, four and five 
sites by presenting the fraction of the sky where the re- 
construction is performed with a given error area. This 
figure also shows a significant improvement of the source 
localization (particularly for the 90% error area) as more 
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FIG. 4: Median error angle vs average detector SNR obtained 
with the unmodeled algorithm for all types of injections and 
different network configurations: LHV, LHVA, LHVJ, LH- 
VAJ. 



sites are used for the reconstruction. The best coordi- 
nate resolution is obtained with the 5-site network and it 
is compatible with the field-of-view of most optical tele- 
scopes. 

The coordinate resolution depends also on the wave- 
form morphology and the polarization content of GW sig- 
nals (see for details Section VI B). If reconstructed with 
the least constrained unmodeled algorithm, the SG waves 
with linear and circular polarization have less accurate 
source localization (see Figure [7]). However, the coor- 
dinate resolution can be significantly improved if recon- 
struction is constrained by the source polarization model. 
In general, the reconstruction improves as more accurate 
models, with fewer free parameters, are used. We ex- 
pect, for example, that the template analysis of waves 
from the coalescence of binary neutron stars and black 
holes [45] should result in a better sky localization than 
for the unmodeled burst search. 

Tables [V]|VIII| summarize the results of the analysis for 
different injection signals, source polarization models and 
networks by showing the fit parameters A and A+B+C, 
which correspond to the median error angle for events 
with high and low SNR, respectively. 




FIG. 5: Median error angle for LHV, LHVA, LHVJ and LH- 
VAJ networks (from top to bottom) as a function of source 
coordinates (0 - latitude, <fi - longitude) for injections with the 
network SNR between 10 and 30. 



B. Waveform Reconstruction 

The signal waveforms are obtained from the solution 
of the likelihood functional. Those are the reconstructed 
detector response as defined by Eq. 3.9 To characterize 
the algorithm performances on waveform reconstruction 
we consider the normalized residual energy A: 



A = 



(£h|£ h ) 



(5.3) 



where the inner products are defined by Equation |2.5 
The same as for the coordinate reconstruction, the accu- 
racy of the waveform reconstruction strongly depends on 
the strength of detected signal and the waveform mor- 
phology. At low SNR, the reconstruction is affected by 
the detector noise, however it improves with the increas- 
ing value of SNR (see Figure |8|. At high SNR, the recon- 
struction accuracy reaches a limit due to the finite preci- 
sion of the algorithm. In this paper we do not present a 
detailed study of the waveform reconstruction. However, 
such work is planned in the future. 
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FIG. 6: Fraction of the sky (vertical axis) for 3- site (LHHV), 
4-site (LHVA and LHVJ) and 5- site (LHVAJ) networks where 
sources are reconstructed with a given 50% (top) and 90% 
(bottom) error region (horizontal axis). 
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WNB HF 
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TABLE V: Pointing accuracy (A+B+C / A) for unmodeled 
search. 



elliptical 


LHV 
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LHVAJ 


SGQ9 LF 


5.3°/0.9° 


1.2°/0.4° 


1.4°/0.5° 


0.9°/0.4° 


SGQ9 HF 
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SGQ3 LF 


3.6°/0.6° 
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0.8°/0.4° 


0.5°/0.4° 


SGCQ9 LF 


8.2°/0.7° 


1.9°/0.4° 


1.5°/0.4° 


0.9°/0.4° 


SGCQ9 HF 


7.2°/0.8° 


1.4°/0.4° 


l.l o /0.4° 


0.9°/0.4° 
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FIG. 7: Median error angle vs SNR for LHV network com- 
paring different constrain searches. Top: SGQ9 LF wave- 
form with unmodeled (black) and linear (red) searches. Bot- 
tom: SGCQ9 LF with unmodeled (black) and circular (red) 
searches. 



linear 


LHV 


LHVA 


LHVJ 


LHVAJ 


SGQ9 LF 


3.670.5° 


1.170.4° 


1.170.4° 


0.770.4° 


SGQ9 HF 


4.5°/0.6° 


0.8°/0.4° 


0.9°/0.4° 
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4.5°/0.5° 


0.8°/0.4° 


0.9°/0.4° 


0.5°/0.4° 



TABLE VII: Pointing accuracy (A+B+C / A) for linear 
search. 



circular 


LHV 


LHVA 


LHVJ 


LHVAJ 


SGCQ9 LF 


2.3°/0.6° 


0.8°/0.4° 


1.070.4° 


0.7°/0.4° 


SGCQ9 HF 


1.7°/0.6° 


0.8°/0.4° 


0.6°/0.4° 


0.6°/0.4° 



TABLE VI: Pointing accuracy (A+B+C / A) for elliptical 
search. 



TABLE VIII: Pointing accuracy (A+B+C / A) for circular 
search. 
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FIG. 8: Normalized residual energy A network SNR ob- 
tained with the unmodeled algorithm and the LHV network 
for all types of injections 



VI. FACTORS LIMITING RECONSTRUCTION 

In this section we describe the factors which limit the 
accuracy of the coordinate reconstruction including: a) 
angular and strain sensitivity of the detectors, b) polar- 
ization content of the signals, c) calibration uncertainties 
and d) limitations of the reconstruction algorithm. 



A. Antenna patterns and detector noise 

The angular and the strain sensitivity of the network 
is fully characterized by its noise-scaled antenna pattern 
vectors (see Eq. 2.1). Due to unfortunate sky location or 
elevated detector noise, the components of these vectors 
corresponding to an individual detector may be small 
with respect to the other detectors, effectively exclud- 
ing this detector from the detection and reconstruction 
of a marginal GW signal. For example, for a source at 
(0 = -40°, <j> = 50°) the angular sensitivity of the VI 

' 0) and it can not 



detector is too low (y/|F + | 2 + \F X 
contribute to the reconstruction unless the GW signal is 
very strong. For the LHV network it means that for a 
significant fraction of the sky the triangulation is per- 
formed with only two detectors, which is not sufficient 
for the accurate source localization. For this reasons it 
is very desirable to have four or more detectors in the 
network operating in coincidence. 



B. Waveforms and polarization 

For a given direction to the source, the reconstruc- 
tion accuracy may be very different depending on the 
signal polarization. For example, a linearly polarized 



signal may not be measured by one of the de- 

tectors for some values of the polarization angle when 
\F+\ is almost null. As a result, with the 3-site networks 
the source localization for such signals is expected to be 
poor (see Figure [9]). On contrary, signals with two polar- 
ization components can be localized well via their cross 
component even at the minimum of F+. 
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FIG. 9: Dependence of (top plot, at — —30° and 
cj) — 72° ) and the reconstruction accuracy on the polariza- 
tion angle ip for two types of waveforms: linearly polarized 
SGQ9 LF waveforms (middle) and randomly polarized WNB 
LF waveforms (bottom). 



C. Calibration uncertainties 

The coordinate reconstruction may be affected by the 
calibration uncertainties of recorded data streams. Typi- 
cally the calibration uncertainties are of the order of 10% 
in the amplitude and few degrees in the phase [40] . These 
systematic distortions of the GW signal may result in a 
systematic shift of the reconstructed sky location away 
from a true source location. 

To estimate this effect we introduce a variation of the 
amplitude and phase of the injected detector responses. 
The amplitude variation is selected randomly between 
±10% or 0% and the random time shifts of ±32 jis or 
fis are introduced. The non-zero time shifts correspond 
to a phase shift of ±2.5° and ±11.5° at the low and high 
frequencies respectively. Such "misscalibration" is ap- 
plied to all detectors in the network. The results are 



reported in table IX for the LHV network. They show 
that depending on the signal morphology and bandwidth 
the calibration uncertainties may affect the coordinate 
reconstruction. The effect of calibration uncertainties is 
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particularly visible at high SNR where the angular resolu- 
tion is less affected by the detector noise. Similar studies 
for the LHVAJ network do not show, even at high SNR, 
any significant impact of the calibration errors used in 
the analysis. This is an expected result, because a better 
constrained LHVAJ network provides more robust source 
localization than the LHV network. 



Waveform 




Amplitude 


Phase 


WNB LF 


0.7° 


1.0° 


0.9° 


WNB HF 


0.4° 


0.6° 


0.8° 


SGQ9 LF 


0.7° 


2.8° 


1.2° 


SGQ9 HF 


0.9° 


1.6° 


1.4° 


SGQ3 LF 


0.5° 


2.1° 


1.0° 


SGQ3 HF 


0.4° 


1.1° 


1.0° 


SGCQ9 LF 


0.8° 


2.5° 


1.1° 


SGCQ9 HF 


1.1° 


1.9° 


2.0° 



TABLE IX: Pointing accuracy (fit parameter A) for the 
LHV network and different signals (column 1): no calibra- 
tion errors (column 2), amplitude and phase mis- calibration 
(columns 3 and 4 respectively). 



D. Reconstruction algorithm 

There are several factors limiting the accuracy of the 
coordinate and waveform reconstruction due to the cWB 
algorithm. For high SNR events the coordinate resolu- 
tion is limited by the cWB sky segmentation which is 
0.4 x 0.4 degrees. Therefore the error angle can not be 
less than 0.4 degrees. Also for the high frequency events 
the coordinate resolution is limited by the discrete time 
delays (see section III) with the step of 1/16384 sec- 
onds and by the accuracy of the time delay filter (few 
percent) used in the analysis. Also in the analysis we 
did not use any unmodeled constraint specific for indi- 
vidual networks, which, in principle, may improve recon- 
struction. These limitations are not fundamental and the 
algorithm performance can be improved in the future. 



is the network response across the sky. The alignment 
factor, which strongly depends on the number of detec- 
tors and the orientation of their arms, determines the 
relative contribution of the two GW polarizations into 
the total network SNR. 

It requires several non-aligned detectors (preferably 
more than three) for a robust detection and reconstruc- 
tion of both GW components. The coordinate recon- 
struction strongly depends on the signal waveforms, net- 
work SNR and the number of detector sites in the net- 
work. The reconstruction can be significantly improved 
when it is constrained by the signal model. Although a 
crude coordinate reconstruction (ring in the sky) is pos- 
sible with the networks of two spatially separated sites, 
at least three detector sites are required to perform the 
source localization. The accuracy of the localization dra- 
matically increases for networks with more than three 
sites, particularly for the low SNR events. For exam- 
ple, the LHHV and LHVA networks are expected to have 
about the same detection rates, however, the 4-site LHVA 
network would have much better performance for the ac- 
curate reconstruction of GW signals. The pointing res- 
olution required for joint observations with the electro- 
magnetic telescopes is achievable with the networks con- 
sisting of four sites. The LHVAJ network demonstrates 
further improvements, both in the detection and recon- 
struction of GW signals, reaching a sub-degree angular 
resolution. In addition, due to the limited duty cycle of 
the detectors, both the LCGT and the Australian detec- 
tors will significantly increase the observation time when 
any of 4-site networks are operational. 

The advanced LIGO and Virgo detectors are very ca- 
pable of the first direct detection of gravitational waves. 
However, for better reconstruction of the GW signals 
more detectors are required. Extra detectors introduce 
an important redundancy which lower the impact of lim- 
ited duty cycle of the detectors, makes the coordinate 
reconstruction more accurate, and less dependent on the 
waveform morphology and calibration uncertainties. The 
construction of the LCGT and the detector in Australia 
will significantly enhance the advanced LIGO- Virgo net- 
work and these detectors will play a vital role in the 
future GW astronomy. 



VII. CONCLUSION 

In the paper we present the results of the source lo- 
calization and reconstruction of GW waveforms with the 
networks of GW interferometers. For a general character- 
ization of the detector networks we introduce few funda- 
mental network parameters, including the effective noise, 
and the network antenna and alignment factors. The ef- 
fective power spectral density of the network noise deter- 
mines the average network SNR for a given population 
of GW signals. For each direction in the sky the network 
performance is characterized by its antenna and align- 
ment factors. The antenna factor describes how uniform 
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